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Exemplo: 


65 carcaças de 
motores electricos 
de um peso bruto de 
500 quilos inteira- 
mente decapados 
em 8 minutos 
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Peso das peças 

por carga 240 Kg 
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à tesoura nada resolve... 


quando se necessita de corrente contínua 


Para tal podemos fornecer-lhe 

outras espécies de «tesouras», graças 
ao nosso vasto programa de recti- 
ficadores. Dispomos de uma fábrica 
de semicondutores, apta a fornecer 
rectificadores para as mais variadas 
aplicações. Os resultados obtidos 
nesta fábrica são cada vez mais 
surpreendentes, devido a constantes 
trabalhos de pesquisa e de aperfei- 
çoamento. Os processos de fabrico 
altamente desenvolvidos e a longa 
experiência adquirida na forma de 
produção mais adequada, garantem 
uma qualidade cada vez melhor. 
Quer se trate de rectificadores de 
selénio com características de 
passagem melhorada, quer dé díodos 
de silício ou de células de silício 
comandáveis (tiristores) com tensão 


de bloqueio aumentado, sempre 

que haja necessidade de corrente 
contínua ou de «moldar», corrente 
alternada, os nossos elementos 
semicondutores e instalações rectifi- 
cadoras constituem a solução. 
Podemos fornecer quase todos os 
tipos de rectificadores normais em 
curto espaço de tempo, do nosso 


armazém, 
medir controlar regular 
automação 
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UMA MANUTENÇÃO DA MAQUINARIA MAIS ECONÔMICA 


Se você evitar por todas as formas ao seu alcance as possiveis falhas da sua máquina, resultará muito 
mais econômica a sua manutenção. 


O óleo industrial ESSO proporciona uma segurança absoluta no rendimento, reduzindo ao minimo 
as avarias. 


O óleo industrial ESSO torna possivel um menor desgaste nas diversas peças, uma melhor conservação 


de todo o tipo de máquinas, o que, afinal, pressupõe um aumento da produção e a obtenção de maiores 
beneficios. 


ESSO... proGRÊSSO.. sucÊSsSO 


O SIMBOLO ESSO REPRESENTA SETENTA E SETE ANOS DE EXPERIÊNCIA 
MUNDIAL, COM VINTE DE BONS SERVIÇOS EM PORTUGAL 


ESSO STANDARD PORTUGUESA (COMPANHIA DE PETROLEOS) S.A. R.L. 
Em Lisboa: Rua Filipe Folque, 2 - No Porto; Av. Serpa Pinto, 227 - Matosinhos 
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Com os condensadores MP os motores de inducao 
bem e funcionam melhor ! 


arrancam 


Os condensadores MP apresentados 
pela ITT em caixas seladas de alu- 
minio são auto-regeneradores. Ofere- 
cem, portanto, uma elevada resis- 
tência aos usuais factores de destruição 
tais como humidade elevada, picos 
de tensão, etc.. 


O alto nível de qualidade tornam-nos 
especialmente indicados para o arran- 
que e funcionamento com motores 
monofásicos de indução. É o caso, 


componentes electrónicos 


por exemplo, dos motores utilizados 
nas máquinas de lavar, queimadores 
de óleo, bombas, ventiladores, serras 
circulares e máquinas-ferramentas 


Características : 


— Capacidade : 0,5 uF a 160 uF; 

— Tensão nominal : 220 V a 600 V; 

— Temperaturas de serviço -—259C 
a 85ºC 


Para informações adicionais relativas 


aos nossos condensadores MP ou à 
nossa linha completa de componentes, 
que se estende desde os micro- 
módulos às válvulas gigantes de 
emissão, consulte a : 


STANDARD ELÉCTRICA, s.a.r.l. 
Avenida da India 

LISBOA 3 

Telefone : 63.81.71 

Telex : 572 


ITT 
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Comunicações 


equipamentos 


e sistemas completos 
+ Cuidadosamente 
projectados 


* Amplamente 
experimentados 


* Completamente 
integrados 
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ASPECTOS DA ANÁLISE DE CORRELAÇÃO 
E DE REGRESSÃO 


RESUMO 


Focam-se aspectos da análise de correlação e de regres- 
são, tratando-se em especial da regressão linear a duas ou 
mais variáveis. Abordam-se problemas de estimação, de 
ensaios de significância e de intervalos de confiança. 


1 —- CONCEITOS DE CORRELAÇÃO E DE 
REGRESSÃO SIMPLES E MÚLTIPLA 


Diz-se que existe correlação entre duas variá- 
veis aleatórias Y e X quando a lei de distribui- 
ção de qualquer delas, para um valor constante 
da outra, depende desse valor. 

O estudo da correlação entre duas variáveis 
(correlação simples) faz-se a partir das distribuições 
bidimensionais. 

Dadas duas variáveis aleatórias continuas Y 
e X, a probabilidade de se ter simultâneamente 
Y<Y<yt+dy ex<X<x+ dx é dada 
pela probabilidade elementar f (y, x) dy dx em 
que f (y,x) é a função densidade de probabilidade 
conjunta do par (y, x). 

A distribuição marginal de Y (ou de X) é a 
distribuição desta variável considerada indepen- 
dentemente da outra. A função fi (y), função 
densidade de probabilidade marginal de Y, é tal 
que f; (y) dy representa a probabilidade ele- 
mentar de ser y< Y <. y + dy, qualquer que 
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por ANTÔNIO QUINTELA 
Eng.º Civil 
Assistente do 1.858,71. 
Bolseiro do 1. A.C. 


SYNOPSIS 


Aspecte of correlation and regression analysis are dealt 
with, particularly those of linear regression involving tao 
or more variables. The paper covers problems of estimation, 


significance tests and confidence limita, 


seja o valor de X. Define-se de igual modo 
fa (x), função densidade de probabilidade mar- 
ginal de X. 
É, pois, 
fi (y) = E (y, x) dx 
fa (x) = / f (y, x) dy 


em que os integrais são estendidos ao domínio 
de variação de X e de Y, respectivamente. 

A distribuição de Y condicionada por X é a 
distribuição da variável Y sob a condição de 
ser x<X< x + dx. 

A função f(y/x), função densidade de proba- 
bilidade de Y condicionada por X, é tal que 


= 


representando f(y/x) dy a probabilidade de 
y<Y< y + dy sob a condição de ser 
XxX < x + dx. 
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Para variáveis independentes, será: 


E (y,x) = fa (x) fi (y) 
É (y) = fi (y) = (y/2) 
f (x) = fa(x) = (x/y) 


As médias e as variâncias das distribuições 
marginais das variáveis Y e X são respectiva- 
mente 


ty==/y ti(y) dy = y f(y,x) dy dx 

Ex = [x fa(x) dx =// x f(y,x) dy dx 
1, =[(y-vy) fi (y) dy = ([(y-uy) E (y,x) dy dx 
02 [(x-pyx)? fa (x) dx = //(x-py)? E (y,x) dy dx 
Definem-se covariância — covyy — e coeficiente 
de correlação — p — entre as variáveis Ye X por 
meio das igualdades 


covyx = // (y-ty) (x-px) F(y,x) dy dx 


— COVyx 


o Gy Tx 
O coeficiente de correlação assim definido é 
—1< (<1 
No caso das variáveis Y e X serem indepen- 
dentes, o coeficiente de correlação é nulo, pois 
de ser f(y,x) = fi (y). fa (x) resulta 


COVyx =0 €0=0 


A recíproca não é verdadeira, pois ;=-o não 
implica 
E (y,x) = fi (y). fa (x) 


O valor médio de Y condicionado por X, ou 
seja, a média dos valores da variável Y para o 
valor de X compreendido entre xe x +dx,e a 

«A - Vc 
variância de Y condicionada por X represen- 
2 —- 
tam-se por H'yjx e 7,,, € definem-se pelas relações 


Hyix = [y f(y/x) dy 


ye sf (y-pyix)? E (y/x) dy 


Quando existir correlação entre as variáveis, O 
valor médio de Y condicionado por X, assim de- 
finido, será função de x. A dependência entre 
Hyix e x é expressa pela equação de regressão de Y 
sobre X e a sua representação no plano (y,x) é 
designada por curva de regressão de Y sobre X. 
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De modo semelhante, se definem pxiy € Srys 
respectivamente, valor médio e variância de X 
condicionados por Y, e também equação e curva 
de regressão de X sobre Y. 

As curvas de regressão de Y sobre X e de X 
sobre Y são em geral diferentes. 

O significado da curva de regressão de Y sobre 
X (Fig. 1) é o seguinte; 


Para um valor de X compreendido no inter- 
valo x; x + dx, os valores de Y apresentam uma 
distribuição (distribuição de Y condicionada 
por X) de média |y;x, situada sobre a curva de 
regressão, e de variância Six que caracteriza a 
dispersão dos valores de Y em torno de Hyjx. 


probabilidade da distribuição 


Curva densidade de 
de Y condicionado por X 


cal Myix= Myuta) 


Fig. 1 


Estes conceitos continuam a aplicar-se no caso 
de variáveis discretas, substituindo-se, nas ex- 
pressões anteriores, os integrais por somatórios. 

Quando se faz referência a correlação, enten- 
de-se que ambas as variáveis Y e X são alea- 
tórias, tornando-se possível a partir do valor de 
uma delas conhecer o valor médio e a dispersão 
dos correspondentes valores da outra, e recipro- 
camente; podem estabelecer-se, portanto, equa- 
ções de regressão de Y sob X e de X sobre Y, 

Pode, porém, estudar-se a regressão de uma das 
variáveis sobre a outra quando esta não for uma 
variável aleatória, não se tratando neste caso de 
correlação entre as duas variáveis. 

Por exemplo, o estudo da dependência entre a 
resistência à compressão de cubos de betão e o 
tempo de cura faz-se a partir do conceito de 
regressão; o tempo de cura, tomando valores 
fixados prêviamente, não pode ser considerado 
como variável aleatória. Já poderá existir corre- 
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lação entre a resistência à compressão e a resis- 
tência à flexão de cubos de betão. 

Os conceitos de correlação e de regressão são 
generalizáveis ao caso de mais de duas variáveis 
(correlação múltipla e regressão múltipla), apresentan- 
do-se a seguir o conceito de regressão múltipla, 

Sejam as variáveis continuas X,, X,, Xas 004 Am, 
tendo a variável Xy uma distribuição condicio- 
nada por X,, X;,..., Xm dependente dos valores 


Xgs X30++ +» Km, com média E ixo x3,..., xm € COM 


xm 
variância 9y,/x,, x3,..., xm" Designa-se por equa- 
ção de regressão de X, sobre X, X,...,X, à 
equação que relaciona L'x,/x,, x;, . 
valores x3, X3, ..., Xm 


-. Xm Com os 


Uxixo,X3, cr ET GO O, 0.04 Xu Db; O ooo DM) 


em que 04, 04, 924,..., x são parâmetros. 

No espaço (X4 Xa,..., Xm), a equação ante- 
rior é representada por uma hiper-superfície que 
se designa por superfície de regressão de X 
sobre Xs, X3,..., Xm (Fig. 2). 


X1 


Superfície de 


Fig. 2 


2 — REGRESSÃO. MODELO E ESTIMAÇÃO 
DE PARÂMETROS 


O modelo estatístico geralmente usado na 
análise de regressão é do tipo seguinte: 


— A média Px/Xy o X xm de X1 con- 


Dada 
dicionada por Xs2, X3,..., Xm é função 


de X,, Xj, -.., Xy é de um certo número 


de parâmetros 1, 02,..., Ox 
*s Xm ne é (Xg, X3 4 ea Xm + 


01, 0a, .. 4 94) 
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— À variável X4 condicionadapor Xsa, Xs,.. 
distribui-se em torno de Px ix, 


«4 Xm 


.0. 4 Xm 
segundo uma lei normal com variância condicio- 
nada constante 


1X» 


2 


o | 2 
XX, 3X voc 0, X od 


Assim, para um grupo de valores das variáveis 
(Xyjs Xojy Xyj, «e» Xmj) à função densidade de pro- 


babilidade de X4 condicionada por X2, X3,...,Xm 
toma o valor 


É (x5;/Xa » X3je ce+s Xmj) = 
e o e 
20º “ai 


2 
— 40% Xyo coco Xmis 84, 0a, 048] | 


1 
— CC C— €X 
vos 


Dada uma amostra de n grupos de valores 
das variáveis (X,, Xs;, Xyjs «--» Xmi), com i= 
=1,2,..., no critério de máxima verosimi- 
lhança para a determinação dos estimadores 01, 
0a, ..., dk dos parâmetros da equação de regres- 
são impõe que seja máximo o produto 


É Callas Ragesey Dn 


: É Du Mapa cs Mailo cos 


E 


1 Y 
ex — —— da XxX = 
Ê 2 q? [a 


co Xmjr Bye Ogrecos do | | 


Esta condição equivale a tornar mínima a 
soma 5 dos quadrados dos desvios entre os va- 
lores de Xi observados e os valores médios con- 
dicionados que se calculam pela equação de 
regressão 


RR. 
(7 V2 nn 


O (Xyjy Xypre 


S=: 


— b (x 1, Xv cos X 


mi” 


b, Duo so dq | 


Os estimadores 01, 13,..., dx são tais que 
minimizam a função S 


dS 
o 04 


d 
Des 
o 0a O 9x 
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O método descrito para a estimação dos pará- 
metros de uma equação de regressão é conhe- 
cido por método dos mínimos quadrados, 

Ao definir o modelo de regressão não se exige 
que X2, X3,..., Xm sejam variáveis aleatórias. 
Apenas se estabelece que a distribuição condi- 
cionada de X4y para valores determinados de 
Xa, X3,..., Xm Seja normal e com variância 
constante. 


3 — DISTRIBUIÇÃO BINORMAL. REGRESSÃO 
SIMPLES LINEAR 


A função densidade da distribuição bidimen- 
sional normal (distribuição binormal) escreve-se 


1 
f x) aci ESSES = —— 
ty 2% Oy 0x1 — q? 
né 2 
E | 2 [Pesa = 
2(1-P)L 9, 


— 22 Guy) (epa) | Se | | 
Oy Ox o, 


A distribuição binormal assim definida é tal 
que: 

— À distribuição marginal de Y é uma distri- 
buição normal com média &y e variância os. 

— A distribuição marginal de X é uma distri- 
buição normal com média x e variância cf. 

— O parâmetro o é o coeficiente de correlação 
entre Ye X. 

— A distribuição de Y condicionada por X é 
uma distribuição normal com média yyjx e variân- 
cia 75, dadas por 

Eylx == ty + E (x-tx) 
Tx 


e RP 
aos Ep) 


— À distribuição de X condicionada por Y é 
uma distribuição normal com média 1x/y e variân- 


cia 7;,, dadas por 


x/y 
Tx 

Exly = Ex + p— (y-2'9) 
Ty 


E ramo 
Daly a 7, (1 pt 22) 


Verifica-se assim que a equação de regressão 
de Y sobre X é linear e que a variância de Y con- 
dicionada por X é independente do valor da va- 

“+ A . - 4 
riável X. Idêntica observação se poderá fazer em 
relação à regressão de X sobre Y. 
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Quando for p ==+ 1, coincidirão, como se 
pode demonstrar, as equações de regressão de Y 
sobre X e de X sobre Y, e as variâncias condi- 
cionadas a e Ty serão nulas, o que significa 
que deixa de se estar em presença de uma rela- 
ção estatística para ter lugar uma relação funcional. 

Quando é :==0, a equação de regressão de 
Y sobre X é representada por uma recta para- 
lela ao eixo dos xx, o que significa que o conhe- 
cimento de x não traz qualquer contribuição 
para o conhecimento de y. No caso particular da 
distribuição binormal, »: = 0 implica a indepen- 
dência das variáveis Y e X. 

Verifica-se a seguir que o modelo de regressão 
de Y sobre X (em que se considere linear a 
equação de regressão, e normal e com variância 
constante a distribuição de Y condicionada por X) 
coincide com o modelo da correspondente re- 
gressão obtida a partir da lei binormal. 

Seja dada uma amostra de n pares de valores 
(yi, xi) das variáveis aleatórias Y e X, perten- 
centes ao universo (y, x) em que a equação de 
regressão de Y sobre X seja linear (uyjx == 
= So + bx) e a distribuição de Y condicionada 
por X seja normal e com variância constante. 

Se se designarem por bo e b as estimativas, 
obtidas a partir da amostra, dos parâmetros fo e 
& da equação de regressão (em que o parâme- 
tro É é designado por coeficiente de regressão), 
esta pode escrever-se 

y — bo +- bx 
em que y' substitui uyix, de que y' é estimativa. 

As estimativas das médias, y e x, e das va- 
riâncias, s; e op das distribuições marginais e 
da covariância, Ea obtidas a partir da amos- 
tra, calculam-se pelas expressões 


E 1 1 e 
pa - o 
y n 74 Sy n 20 y) 
ad a EM = 
1 mi 
ga — D) (yi— y) (x; — X) 
y n | 


O coeficiente de correlação é estimado por 
a (yi—y) (xi — x) 


p=: 
n Sy Sx 
As estimativas bo e b dos parâmetros £, e £ 


podem ser obtidas, nas condições enunciadas, 
pelo método dos mínimos quadrados 
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bo <a Yo bx 
2 (yi—y) (xi—x) 
A ii at 
2 (x,—x)º Sx 
à 
escrevendo-se a equação de regressão 
A = 
y=y+r(x—s) 
Sx 
A estimativa da variância de Y condicionada por 
X, suposta por hipótese independente do valor 
de X, obtém-se a partir dos desvios entre os 
valores observados y, e as correspondentes esti- 


mativas dos valores médios condicionados de Y, 
situados sobre a equação de regressão, y*, 


& (y; — y'9* 


n 


DE a, que 
né = Syjx 


e demonstra-se, tendo em conta as expressões 
anteriores, que é igual a 


s=s iu =s, (1 — 1º) 


yr 

Observa-se, assim, a concordância das ex- 
pressões agora obtidas para y' e Sox com as obti- 
das a partir da distribuição binormal. 

Para avaliar o significado da estimativa do 
coeficiente de correlação, veja-se que a expressão 
do seu quadrado 

| Z(yi—y) u—x) |? 
3 O Aa 


n Sy Sx | 


se transforma, no caso de regressão linear, nas 
expressões seguintes : 


Zy—-y)?—-Z(y—y)? 
É RE. SAE fo aa he 


: (yi — y)? 


z (y;— y)? 


2 (y—y? 


Da última expressão sobressai o significado 
do coeficiente de correlação: o seu quadrado 
indica a fracção da variância de Y que é expli- 
cada por X, através da equação de regressão 
linear. 

Deve notar-se que a estimativa r do coefi- 
ciente de correlação (ou o seu quadrado) dá uma 
indicação da intensidade da dependência entre 
Y e X sômente na hipótese da regressão linear. 
Pode obter-se um valor de r vizinho de zero, 
sendo os valores de Y explicados por X através 
duma equação de regressão doutro tipo. No caso 
da Fig. 3, r toma um Valor próximo de zero 
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Fig. 3 


em virtude de terem sinais contrários os pro. 
dutos (yi—y) (xi—x) correspondentes a pon- 
tos situados acima (ou abaixo) da ordenada y 
e a um e outro lado da abcissa x. 

A estimativa s? de variância de Y condicionada 
por X calculada a partir de uma amostra de n 
pares de valores, deverá ser corrigida para for- 
necer uma estimativa não enviesada, s Ml da 
variância condicionada de y do universo a que 
pertence a amostra. Para obter s””, a soma dos 
quadrados dos desvios y;— y'i deverá ser divi- 
dida, não pelo número de desvios, mas pelo nú- 
ro de graus de liberdade, ou seja, pelo número 
de desvios independentes. Como a determinação 
de y'; necessita de dois parâmetros da equação 
de regressão, o número de graus de liberdade 
é de n— 2 

| a (yi— y'9)? 


n 
Ss = ———— = — 


n—2 n— 2 


sy (1 = 1º). 


Também a estimativa corrigida do coeficiente 
de correlação será tal que 
=D 


n—2 


As expressões e conceitos apresentados para o 
modelo da regressão linear de Y sobre X conti- 
nuam a ter validade ainda que X não seja uma 
variável aleatória. Simplesmente, se X não for 
variável aleatória, os parâmetros x e sx serão a 
média aritmética e o desvio quadrático médio das 
observações xi; e não as estimativas da média e 
do desvio-padrão da distribuição marginal de X; 
também Sox não será a estimativa da cova- 
riância, 
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4 — REGRESSÃO SIMPLES LINEAR. DISTRI- 
BUIÇÕES DAS ESTIMATIVAS E INTERVA- 
LOS DE CONFIANÇA 


No número anterior, tratou-se da regressão li- 
near e calcularam-se, a partir de uma amostra de 
n pares de valores, as estimativas dos parâmetros 
da equação de regressão e do coeficiente de corre- 
lação. 

Importa investigar em que medida os verda- 
deiros valores dos parâmetros poderão diferir das 
suas estimativas pela acção das flutuações devi- 
das ao acaso. Em particular, interessa investigar 
a probabilidade de as estimativas b e r do coefi- 
ciente de regressão É e do coeficiente de corre- 
lação » assumirem os valores calculados, quando 
a amostra fosse proveniente de um universo (v,x) 
em que $=o g== 0. Interessa ainda determinar 
os intervalos de confiança das estimativas calcu- 
ladas por meio da equação de regressão. 

A estimativa b do cueficiente de regressão tem 
a variância estimada por 


Por outro lado, a variável 


-b-8 


Gb 


t 


segue a distribuição t de STUDENT, com n-2 graus 


de liberdade 


É possível, portanto, determinar o intervalo 
que compreenderá a um nível de probabilidade « 
o verdadeiro valor de £ 


b “+. ta 9 Sb 


em que z é a probabilidade de a variável t ex- 
ceder em módulo tx/a. 

É ainda possível ensaiar a hipótese de $=o 
para determinar com que nível de probabilidade 
x se poderia obter, na amostra proveniente de 
um universo em que fosse É == o, uma estimativa 
do coeficiente de regressão maior ou igual a , 
por flutuações da amostragem. Basta notar que 
a variável t se reduz neste caso a 


| T Vn-2 


t= — 


Sb V 1-rº 
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determinando-se, pela tabela da distribuição de 
SrupenT, a probabilidade x com que o valor cal- 
culado de t é excedido, em módulo, com n-2 
graus de liberdade. 

Se « for pequeno (normalmente estabelece-se 
um valor limite de 0,05), rejeita-se a hipótese 
4-==o. O valor de « é o nível de significância com 
que o coeficiente de regressão difere de zero. 

Em vez de adoptar o procedimento descrito é, 
por vezes, mais cómodo consultar uma tabela 
(MoRice & CHARTIER, p. 565) que indique direc- 
tamente valores de r correspondentes a diferen- 
tes probabilidades de t, e que foi construída por 
FismeR de acordo com a expressão anterior. 

Apresentam-se a seguir as distribuições do 
coeficiente de correlação. 

Quando o coeficiente de correlação do universo 
é nulo (p =), o valor 


o rv n-2 


VI1-rê 


segue a lei de STUDENT com n-2 graus de liber- 
dade. 

Comparando esta expressão com a anterior, 
conclui-se que são equivalentes os ensaios das 
hipóteses $=0 e0=0, 

Para a determinação do intervalo de confiança 
em que deverá situar-se o coeficiente de correla- 
ção do universo p, terá de recorrer-se à variável 
transformada de FisHER 


que segue uma lei sensivelmente normal de mé- 
dia vz e variância 7, dadas por 


ur 


ie a 
E a MY Log 


A média v, da variável transformada estará 
compreendida, a um nível de probabilidade «, 
no intervalo 

za U,,/2 *; 
em que u,, é o valor da variável normal redu- 
zida com probabilidade 2 de ser excedido em 
módulo. 
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Dos limites do intervalo de confiança de p; 
obtêm-se os correspondentes limites do coefi- 
ciente de correlação. 

A transformação de FisnkR pode ser ainda 
utilizada para investigar se é significativa a dife- 
rença das estimativas r, e r, de coeficientes de 
correlação, obtidas de amostras com n, e nº 
pares de valores. Para isso ensaia-se a hipótese 
o ==p2 (coeficientes de correlação dos universos 
iguais). A diferença z, — z, seguirá uma lei nor- 
mal com média nula e variância 


1 + 1 


n— 3 n; — 3 


Calcula-se, portanto, o valor da variável nor- 
mal reduzida 


Quando se ensaia a hipótese py ==23 contra a 
hipótese alternativa p1=Fp2, determina-se, por 
meio de tabelas da lei normal, a probabilidade a 
desta variável exceder em módulo o valor 
de u calculado, ou seja, a probabilidade de 
ser ju| > u,,,- 

Quando a hipótese alternativa de py ==p2 for 
oq >e2, (obteve-se, portanto, r, =>r,), determi- 
na-se a probabilidade x de a variável normal 
reduzida exceder o valor de u calculado (proba- 
bilidade de ser u>u,). 

Para valores pequenos de «x, rejeita-se a hipó- 
tese estabelecida e aceita-se a alternativa. 

A estimativa y' de »yx obtida pela equação 
de regresão y' ==bo + bx tem a variância esti- 
mada por 


eo 


É (xi — x)?- 


Ássim, para um certo Valor x, o verdadeiro 
de tyx (correspondente à equação de regressão 
do universo Lyx==[o-tHfx) estará compreen- 
dida, a um nível de probabilidade =x, dentro do 
intervalo de confiança 

bo + bx “+ ta)» Sy' 
em que t,;, é o valor da variável de STUDENT 
com probabilidade « de ser excedido em módulo, 
com n—2 graus de liberdade. 
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A variância de um valor individual y, sendo 
a soma das variâncias do valor médio condicio- 
nado y' e da variância condicionada, estima-se 
por 


8 gs” E E A + Em o | 
Po EGg—x? 


Assim, quando a variável X tomar o valor x, 
um valor y da variável Y estará compreendido, a 
um nível de probabilidade «x, no intervalo 


bo + bx e ta/, Sy—y' 


Numa representação gráfica, os intervalos de 
confiança correspondem a faixas para as quais 
aumenta a diferença de ordenadas quando x se 
afasta do seu valor médio x. 

No Quadro 1 reproduz-se um modelo de fo- 
lha de cálculo para a análise de regressão sim- 
ples linear. 


5 — REGRESSÃO MÚLTIPLA LINEAR 


Suponhamos que a distribuição de Xi: condi- 
cionada por Xs, X3,... Xm é normal e que: 


— a média é dada por uma equação de regres- 
são linear 


Exif X0 000, Xm = PL Poxa + 
+ Ba x3 + “+ Bm Xm 


— € à variância é constante: 


Mala 1X 4º ++ Xm 
Os coeficientes 22, &3,..., Om denominam-se 
coeficientes de regressão. 
A partir de uma amostra de n grupos de va- 
lores das variáveis será obtida a equação de 
regressão 


x1= by + baxa + b;xy +... + bm Xm 


em que x'4, by, bs, bs,..., bm são estimativas de 
Exa,» Ryrcc0s Xmr Êa, Êo, o Em. 


No caso da distribuição condicionada de X« 
ser normal, e de a respectiva variância ser inde- 


pendente dos valores das variáveis X2, X3,...1 
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QUADRO 1 


Regressão linear simples 


os -á — valor observado de ni 
y =y+ b(x—x) | y se 
ai x — valor observado de 
y x 
o I — Valores auxiliares 
iy= SZx= 
Ly = EI = o YA = 
1 as 1 
y =: cm E y qu = — py) == 
n n 
(-y)* = (2 x)? = 
=: (y)* 
: n né 
Y id 
, o X |O 
sá a 2%? e E 
n n 
U = 2y. =x 
Sox cs 2yx Es “ nai 
n nº 
2 — Equação de regressão 
2 
b pad “yx — y T— X 
q“ 
" E ps 
3 — Estimativas do coeficiente de correlação r e r' e do desvio-padrão condicionado s' 
+ a 
5 1) rº—1 . EE a 
PP =— — = r? — (n-1) rº—1 G* = a ns,= 
Ss, Sy À —— 2 n— 2 
r= r = S == 


4 — Signiflcância da estimativa do coeficiente de correlação r (x — nível de significância com que r difere 
de zero). 
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Xm, O critério da máxima verosimilhança conduz 
à estimação dos parâmetros pelo método dos mí- 
nimos quadrados, os quais são, portanto, determi- 
nados pela condição de ser mínima a soma 


Z(xi—bi—ba x2— byxs —... — bmXm)? 


em que se omitiu o índice correspondente a cada 
grupo de valores das variáveis. 
Esta condição conduz a 


X =x + ba (x2— x3) + 
+ bs(xs— x3) +... + bm (Xm— xm) 
em que X1, X3, X3,..., Xm São as médias aritmé- 
ticas dos valores das variáveis, e a um sistema 
de m-1 equações com m-1 incógnitas que são as 
estimativas dos coeficientes de regressão 


Os MLS E o RE 
E RE + bm 2 X Xm == & Xi Xa 


b2X, X,+b, EX + 
+ ...+Hbma 2 X3 Xm=2 XX 


ba Z X2Xm + bs 2 X3 Xm + 
+..+bn2 Xi = 2X Xm 


No sistema anterior e nas expressões apre- 
sentadas a seguir, as maiúsculas referem-se a 
desvios das variáveis em relação às suas médias 
(Xi=x—x;). 

O coeficiente de correlação múltipla R é esti- 
mado por 


Eus)! Gu —x)? Elba)? 


Rº= 
(xi — x1)? 3 (x — x1)? 

em que x1 e x representam valores respectiva- 
mente observados e calculados pela equação de 
regressão. 

O coeficiente de correlação múltipla é tal que 
o seu quadrado indica a fracção da variância da 
variável Xi que é explicada por X3, X3,...,Xm 
através da equação de regressão linear. 

Demonstra-se que R? se pode escrever sob a 
forma 
= 
EXi 

+ bs 2X X +... + bm 2X Xm) 


Rº = (br ZXi X3 + 
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A estimativa corrigida R' do coeficiente de 
correlação do universo, de que foi extraída a 
amostra, obtém-se por 


ni (n—1) RR—-(m—1) 
n—m 


R'? 


A estimativa da variância de Xi condicionada 


por X1, X3,..., Xm, tendo em conta o número 
de graus de liberdade, é 
“PR? 
GS! — 1—-R s x 
n—m 


As variâncias das estimativas bs, bs,..., bm 


dos coeficientes de regressão (3, É3,..., Em são 
obtidas por 

E a O do sê 19 à aÃ 12 
Sb, Co 9 Sh, = gg 5 + Sh E Cm? 


A determinação de cada um dos coeficientes 
Cor» C33» ---y Com eXige a resolução de um sis- 
tema de m—1 equações a m — 1 incógnitas. 

Por exemplo, c,, é obtido por meio do sistema 


cos 2X5 + Cy 2X X3 +... +Hcom 2X Xm=1 
Co 2X2Xi + co 2X) +... + Cm 2X3Xm =0 


Co 2X, Xm + Cy 2X3Xm +... + co 2X =0 
a partir do qual se obtêm ainda c,;, ..., CGm- 
Um sistema análogo seria obtido para as incóg- 
nitas Cj, «.. + Cym em que os valores do segundo 
membro das equacões seriam, pela mesma or- 
dem, 0, 1, 0,...,0. O último sistema teria 
como incógnitas Cm2, Cm3, -.., Cmm € OS valores 
dos segundos membros das equações seriam O, 
0,...,1. A resolução destes sistemas é facili- 
tada pela circunstância de se verificar cij = cj. 
As variáveis 


ba — fia 


dis ias 
t=-——— E PI Es bg O ps Bm — Em 
Sb, Sb, Sbm 


seguem a lei t de STUDENT com n —m graus de 
liberdade. 

Podem assim determinar-se os intervalos de 
confiança para os coeficientes de regressão, a um 
nível de probabilidade «a, 


ba + sp, ts), ba + sb, ta/, bm É Sbm tx/> 


em que « é a probabilidade de ser |t| > ty/- 
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Poderá ainda utilizar-se a distribuição t de 
STUDENT para avaliar se os coeficientes de re- 
gressão diferem ou não significativamente de 
zero. Admite-se a hipótese [;=0 (i=2,3,...,m), 
calcula-se 


e aceita-se ou rejeita-se a hipótese consoante 
for elevado ou pequeno o valor « da probabili- 
dade da variável t exceder em módulo o valor 
calculado (n —-m graus de liberdade). 

A significância global da regressão pode ser 
ensaiada atendendo a que, no caso do coeficiente 
de correlação múltipla do universo ser nulo, o 
valor da variável 


segue a lei de Srubext com n—m graus de 
liberdade. Os valores de R correspondentes a 
diferentes probabilidades de t dados pela ex- 
pressão anterior constam da tabela construída 
por FisHEkr (MoRicE & CHARTIER, p. 565). 

A estimativa x, de px; /x,, x, **- am, Obtida pela 
equação de regressão, tem uma variância esti- 
mada 

S'4 - o 
8d geem + Shy (x2— x)? + sp; (xa — x9)2 +... + 
Ê 


9 is 
Ki 9 m(Xm — Xm)º 


Quando as variáveis X2, X3,..., Xm tomarem 
os valores x2, x3,...,Xm, O verdadeiro valor 
Exx,xy---sxm estará compreendido a um nível 
de probabilidade « no intervalo de confiança 


X, Em tz/2 Sx; 
em que t segue a distribuição de SrupexT com 
n-m graus de liberdade. 


A variância de um valor individual xi de Xi é 
estimada por 


, , Sa ' 52 
Sw, dá = 52 E day: “e JA (xa — xa)? r 
+ si (x3—x3) +... FS (xm — Xm)? 
J 
e o intervalo de confiança de um valor indivi- 


dual x1, quando as variáveis X2, X3,...,Xm to- 
marem valores x2, X3,..., Xm, Será 


X| ER a/s Sal =, 
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No Quadro 2 reproduz-se o modelo de folha 
de cálculo para a análise de regressões lineares 
a três variáveis. 

Na escolha de variáveis a introduzir no se- 
gundo membro de uma equação de regressão de- 
verá atender-se às considerações seguintes. 

Pode investigar-se se a introdução de uma va- 
riável suplementar Xm-+1 melhora ou não o ajus- 
tamento (ou seja, se reduz ou não a dispersão 
entre os valores de Xi observados e calculados 
pela equação de regressão) calculando o coefici- 
ente de correlação múltipla (ou o desvio-padrão 
condicionado) para as m-+1 variáveis. Aquela me- 
lhoria tem lugar se o valor do coeficiente de cor- 
relação múltipla para m+1 variáveis for superior 
ao valor do mesmo coeficiente para m variáveis 
(ou o desvio-padrão condicionado for menor), 
Em vez de proceder àquele cálculo, pode reali- 
zar-se a mesma investigação utilizando um pro- 
cesso gráfico, menos trabalhoso. Para isso, mar- 
cam-se no gráfico os pontos correspondentes aos 
desvios 4;=xu — x'j entre os valores da variá- 
vel Xr observados (xi) e calculados (x') pela 
equação de regressão 


x1=> + ba (x3— x2) + 
+ ba (x3 — 3) +...+ bm fine ii) 


e aos valores observados da variável suplementar 
X(m + 1)i- 

A representação gráfica dos pontos permite 
avaliar se os desvios podem ou não ser explicados 
pela variável suplementar. Em caso afirmativo 
(Fig. 4) proceder-se-ia à análise de regressão a 
m+1 variáveis. 

A verificação de correlação intensa entre as 
variáveis do segundo membro da equação de re- 
gressão é inconveniente. Com efeito, se for elevado 
o coeficiente de correlação múltipla Rs.3,,,...,m 
quando se considera a regressão de X> sobre 
X3, Xi... Xm, O erro da estimativa bz do coefi- 
ciente de regressão será também elevado, como 
se verifica se se tiver em conta que a estimativa 
da variância de bz (EztkiEL, p. 322; FORD, p. 21) 
pode ser calculada por 


S'2 
Vwy2 2 
am 3(1—R 2.3, 6, v» 


2 
am e 
E “4 m) 


Assim, não obstante a verificação de um bom 
ajustamento (R elevado), podem obter-se valo- 
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QUADRO 2 


Regressão linear a três variáveis 


= o x1 == valor observado de 
x1==x1 + ba (x2 — x2) + ba (x3— x3) xa == valor observado de 


x3 == valor observado de | 


X1 x2 x3 1 — Valores auxiliares 

2xi= Exj= 2x1 X%34 == 

2x2 = 2x; = 2x x = 

2x3 = 2x2 — 2x X3 == 

-— 2x — Zx2 = 2x3 
X — = x = = 4=— = 

n n n 

(> 4)? == 2 x4 2 x3= 

(2 x2)? = Zx.2Zx= 

(- x3)? ii 2x 2 x3= 


1 e 


ig=24-1 Emt= 


; 1 ' 
sy? — Dia E ES ii 
2X X=2%x—— 2x. x = 
| 1 
2 XX = 2 xx — — da. dx = 


1 - , 
2X X=2xx3— — Exa. 2x3 == 
| A=3X2.5X2-— (E Xa Xo)? = 
2 — Equação do regressão 


| EX Xa.5X — EX1Xs. 2X X 
aii CS 


tro ZM EM E 
A 


3 — Estimativas do coeficiente de correlação R e R' e do desvio-padrão condicionado S' 


y 2 
Rt = Dr2XAXrtbal XX Rei (n-1) R-2.. EA E em 
3 Xi = 3 a 
R = R'º = S' = 
4 — Significância das estimativas dos parâmetros (x — Nível de significância com que as estimativas diferem 
de zero) 
= x — 
2% je dh À 
Ca — eo CM = b=— =] = Ps 
» A 2? S'Vcas 
Y 2 
a Xa V Cs; — b3 
= ———— == Cãa == to. === -—— == 
Cas == 3 b 
e A *  SVes 
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db =n4-*4 


mM =mebo(xg-X9) eba (xg-X3) +... +Dm(Xm -Xm ) 


Fig. 4 


res das estimativas dos coeficientes de regressão 
com sinais contrários aos que faria prever o 
conhecimento do mecanismo físico dos fenóme- 
nos a que se referem as variáveis, como se 
exemplifica a seguir. 


Usualmente, é elevado o coeficiente de corre- 
lação entre escoamentos anuais numa dada sec- 
ção de um rio e as precipitações ponderadas 
sobre a respectiva bacia. 

Para a ribeira de Arade em Casa Queimada 
(Quadro 3) a estimativa daquele coeficiente é 
0,840 e a estimativa b'==0,579 do coeficiente 
de regressão [” dos escoamentos anuais (X1) 
sobre as precipitações anuais (X») difere signi- 
ficativamente de zero (nível de significância 
a < 0,001). 

Admite-se que a concentração da precipitação 
anual em meses consecutivos possa explicar em 
parte os desvios entre valores de X1 observados 
e calculados pela equação de regressão anterior 


x1=b, + b'x: 
Assim, para aquele curso de água calculou-se 


a equação de regressão dos escoamentos anuais 
(X1) sobre as precipitações anuais (X:) e sobre 


QUADRO 3 


Equação de regressão com as variáveis do segundo membro correlacionadas 


Variáveis 

X1 — Escoamento anual expresso em altura de 
água (mm). 

Xs— Precipitação anual ponderada sobre a 
bacia (mm). 

X3 — Precipitação no semestre chuvoso pon- 
derada sobre a bacia (mm). 


Equações de regressão 


(1) Uxilxo nm Eo + É" xa. 

(2) Exi/xs aa Po + É"'xa. 

(3) Bxylxy, xy == By + Baxa + Es xa. 

Estimativas dos parâmetros das equações de re- 


gressão 


(1) bo =-—184 b'+sp' = 
(2) bo==211 b'+spi= 


0,579 + 0,100 
0,752 + 0,070 


Notas : 


Número de observações n = 16 


(3) bi =—186 b2+sp=—0,2974-0,166 
bs+sb;y== 1,075 + 0,192 


Estimativas dos coeficientes de correlação e dos 
desvios-padrão condicionadcs. 


(1) r =0,840 r == 0,828 s'= 84 
(2) r = 0,945 r' == 0,941 s == 51 
(3) R=0,956 R'=0,949 SS =47 


Níveis de significância « com que os coeficientes de 
regressão diferem de zero. 


b' 

(1) t= — =5,8 a< 0,001 
Sb 
b” 

(2) t= — = 10,8 «< 0,001 
Sb” 

(3) t= UE qm 1,8 «=0,10 
Sb? 

bs 
= —— = 5,6 «a < 0,001 

Sby 


Estimativa do coeficiente de correlação entre X; e X3: r= 0,941. 
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— 


ás precipitações nos semestres mais chuvosos 
dos vários anos (X3) 


x4 = bi+4- ba x2 + byx3 


A estimativa do coeficiente de correlação é 
mais elevada do que a do coeficiente de corre- 
lação entre precipitações anuais e escoamentos 
anuais. A estimativa bz == — 0,297 do coeficiente 
de regressão [:2, sendo negativa (o que é ilógico), 
não difere significativamente de zero em virtude 
de ser elevado, pela forte correlação entre X3 e Xa, 
o valor da estimativa do seu desvio-padrão con- 
dicionado sb, = 0,166. 

Note-se que para a equação de regressão dos 
escoamentos anuais sobre as precipitações nos 
semestres mais chuvosos dos vários anos 


x1=b/+b”xs 


é positiva a estimativa b” do coeficiente de 
regressão. 
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Quando algumas das variáveis do segundo 
membro da equação de regressão se apresentem 
fortemente correlacionadas, é recomendável com- 
biná-las numa única variável a introduzir na- 
quela equação. 
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CG. D. O. 624.344.5/91 


ELEMENTOS SOBRE A PRODUÇÃO E O CONSUMO DE ENERGIA 
NA REDE ELÉCTRICA NACIONAL **) 


I —Breve nota mensal 


Verificaram-se descarregamentos, sobretudo nos 
aproveitamentos a fio de água e na primeira quinzena 
do mês. 

As afluências ao conjunto do sistema, verificadas 
durante o mês de Maio, têm uma probabilidade de se- 
rem excedidas da ordem dos 59º/,. O coeficiente de 
produtibilidade hidroeléctrica do mês foi de 0,92. 


[1 — Elementos gerais (GWh) 


Acumulados desde 1 de Janeiro 


Variação 
“o 


Produção hidráulica (Ph) ... |1706,2| 2200,5| + 29 
Produção térmica (Pr)... .. | 147,8 O|— 100 
Produção total (PT). ..... 1353,5 | 2200,5/ + 19 
Energia recebida de empresas 


| a) Mensais 
| 1965 | 1966 


não pertencentes ao RNC (Er) | 144 20,5 |+ 42 
Exportações (Ex) .. ...... 1,2 a] — 100 
Importações (1) ......... 97,0 0,0 | — 100 
Saldo importador (Si). ..... +95,5 0,0) — 100 
Consumo em bombagem (Cp)(!) 6,8 46,7 + 587 


Produção para con- é 
sumos perman. (Pcp) (2).. |1736,0| 17964 + 3,5 
Produção para con- É 
sumos não perman. (Pcnp)... | 220, Mdi+ 1 
TOlGLs snacas esse cet 1x 1963,7 | 2221,0 | + 13 


Coeficiente de hidrauricidade | 0,76 1,59 


NOTAS 


(1) Energia consumida em bombagem para armazenamento na Al- 
bufeira do Alto Rabagão, 
(2) A produção para consumos permanentes (Pcp) é determinada 


pela seguinte expressão ; Pep = PrT+4 E, + Si— Ch — Penp 


HI — Diagramas de carga dos dias característicos 


4.* teira: 


1965 | 1966 
Produção hidráulica (Ph) MWh| 12181 11453 
Produção térmica (Pr) MWh 0 Ú 
Produção total (PT) MWh| 12181 lt 53 
Trocas com | Export, (Ex) MWh 0 O 
"Espanha | Import. (1) Mwh| 305 0 
“Uonsum em bomb. hidroel, (Cb) MWh O Ú 


Prod. para cons. perm. (Pp) MWh| 10835 12041 
Prod. para cons. não perm, (Penp) NWh| 1651 24:92 


— TOTAL PTH (|-Ex) MWh| 12486 | 14453 
Potência máx.  MW| 7.06 801 
é á w.,| Potência min,  MW| 320 388 
SE PTH CÃO  Ueiliz, da ponta horas | 17,6 | 18,0 
É E : Factor de carga 0,73 | 0,75 
vê Potência máx. MW 635 | 699 
. p Potência min. MW. 250 265 
Um ar Utiliz. da ponta horas Ni | 17,2 
Factor de carga 0,71 0,72 
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MAIO 


H1 — Diagramas de carga dos dias característicos 


BBZS=« 


7 
; 
g 
- 
= 
o 


e 


ue 
E 


IV — Energia armazenada nas principais albufeiras 


No fim do mês 
Albufeiras : nt q 
GWh % (1) 


Alto Rabagão +... +... | 967,6 99,4 
Paradolês » are us use wc] EMA 100,0 
Venda Nova . « . «o «| 1248 971 
Salamonde , . «+. .... 27,5 99,6 
Caniçada «cocesecsoe 31,9 96,4 
MRE: q do O OIE E E DOS E É 91,1 98,4 
RB «a a suas as co e é) BSS 99,9 
Castelo do Bode. +... ,.| 1628 99,9 
QUO «sura mes 7,3 88,0 
Lagoa Comprida . ..... 211(?) 60,2 
Santa Luzia é c cor css 60,0 97.4 
Pracana . «cc... .. 11,1 86,0 
POOR sa sim es o 12,0(3) 89,4 
Total com A, Rabagão + «| 2077,9 48,5 

sem A. Rabagão . . .| 1110,3 97,8 


NOTAS 

(1) Coeficiente de enchimento em energia das albufeiras, 

(2) Inclui 33 GWh armazenados em Vale do Rossim no início do 
mês e 3,4 GWh no fim do mês. 

(3) Inclui 2,7 GWh armazenados no açude do Poio no início do mês 
e 2,7 GWh no fim do mês. 


(*) Liementos extraídos das estatisticas mensais do Repar- 
tidor Nucional de Cargas (R.N. 0.) As produções e 
os consumos das empresas do R. N. O. representam 
cerca de 94 º/, dos totais do Pais, 
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Meterial fabricado nas nossas oficinas, pronto para expedição 


Fabricação de aparelhagem eléctrica de alta e baixa tensão 
Equipamentos electromecánicos completos para centrais e fábricas 
Postos de transformação tipo monobloco 

Quadros blindados, capsulados e em armário 

firmaduras de iluminação de todos os tipos 

Contadores eléctricos e instrumentos de medida 

Telemedida e telecomando 

Instalações elevatórias de água 

Instalações de saneamento e esgotos 

Ferramentas de corte 

Motores Diesel industriais e maritimos 

Fogões eléctricos 

Frigoríficos 

Ventilação e condicionamento de ar 

Aparelhagem de Raios X e electromedicina TECNICA XY 


ESTAMPARIAS TINTURARIAS 
LAVANDARIAS 


PARA OS VOSSOS TECIDOS TÉCNICOS: 
TRANSPORTADORES, MANCHONS, LAPPING, WOLFRIES, 
FLANELAS, CANEVAS, ETC,, 


EM FIBRAS NATURAIS OU SINTÉTICAS 
CONSULTEM: 


FANAF EL 
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SOCIEDADE LIMITADA 


APART.: 9 
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SONDAGENS geológicas e geotécnicas — CAPTAÇÕES de águas 

subterrâneas — REBAIXAMENTOS do nível freático — CONSOLIDA- 

ÇÕES E ESTABILIZAÇÃO do solo — INJECÇÕES DE CIMENTO 

e oulros produtos — INFRAESIRUTURAS de barragens, luneis e pontes 
e FUNDAÇÕES de todos os tipos 


Única Empresa da especialidade com sede no Ultramar Português 


DELEGAÇÃO - LISBOA SEDE - LUANDA MOÇAMBIQUE 
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TECNICA XVI 


C. D. U. 554.354 (469) 


CONTRIBUIÇÃO PARA O ESTUDO DO ESTUÁRIO DO TEJO 
— SEDIMENTOLOGIA — 


SUMÁRIO 


Com o fim de definir as origens e condições de deposição 
dos sedimentos do estuário do Tejo e da costa atlântica, 
nas suas proximidades, colheram-se cerca de trezentas amos- 
tras distribuidas pelo estuário e praias adjacentes. 

Com base em elementos obtidos de várias fontes, opre- 
sentam-se, em primeiro lugar, considerações relativas à 
distribuição e evolução dos fundos do estuário, à salini- 
dade, ao transporte dos sedimentos e à geologia das forma- 
ções que marginam o estuário e das que constituem a costa 
atlântica adjacente, 

Descrevem-se resumidamente os metodus de estudo apli- 
cados neste trabalho (minerais pesados, carbonatos e gra- 
nulometria) e discute-se a precisão do doseamento dos 
minerais pesados, 

Em consequência dos resultados experimentais obtidos, 
consideram-se : sete associações de minerais pesados, sendo 
três relativas aos sedimentos do estuário e as restantes 
aos sedimentos das praias adjacentes; três grupos de sedi- 
mentos relativamente às percentagens de CaCO,; dois tipos 
granulometricos para os sedimentos das praias allânticas 
e dois tipos granulométricos fundamentais, um arenoso, 
outro sito-argiloso, para os sedimentos do estuário, 

Em face de todos estes elementos, que permitiram definir 
as origens dos sedimentos do estnário e da costa atlântica 
adjacente, apresenta-se una interpretação relativa à evo- 
lução da zona considerada e ao trajecto dos respectivos 
sedimentos, a partir do início da última glaciação do 
Pleistocênico, 


1— INTRODUÇÃO 


Em fins de 1960 o autor teve pela primeira 
vez a oportunidade de trabalhar no domínio da 
Petrografia sedimentar, ao ser-lhe solicitada cola- 
boração num estudo sobre o rio Arade, tendo, a 
a partir de então, vindo a tomar contacto com 
com as técnicas aplicadas ao estudo dos sedi- 
mentos e com a bibliografia da especialidade. 

Dado o grande número de técnicas aplicáveis 
ao estudo dos sedimentos, um dos problemas 
que surgiu, quando do planeamento do presente 
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SYNOPSIS 


In order to define the sources and the depositional con- 
ditions of the sediments in the Tagus estuary and in the 
neighbouring Atlantis const, about three hundred samples 
were collected from the estuary and the adjoining beaches: 

From data obtained from different sources, the author 
discusses, firat, the distribution of and the changes in the 
river bed sediments, the salinity, the transport of sediments, 
as well as the geologi: characteristics of the formation bor- 
dering the estuary and the adjoíning Atlantic coast. 

The author briefly describes the applied methods (heavy 
minerals, carbonat content and grain-size analyses) and 
discusses the aceuracy of the heavy fraction analysis, 

From the experimental resulta obtained, the author con- 
siders: seven heavy mineral associations, three of which 
concerning the estuary sediments and the others the sedi- 
ments of the neighbouring beaches; three groups aceording 
to carbonate content; two grain-size types for the Atlan- 
tic beach sediments and two prêncipal grain-size types, the 
one a sandy, the other a silt-clay type, for the estuary sedi- 
menta, 

From all these data, which made possible the determina- 
tian of the provenance of the sedimenta of the estuary and 
the adjoining Atlantic coast, the author presents an inter- 
pretation of the evolution of the area in study and of the 
movement of the sediments after the early period of the last 
Pleistocene glaciation. 


trabalho, foi o da escolha daquelas que, sendo 
de mais simples aplicação, fornecessem mais úteis 
resultados, tendo em vista os seguintes objec- 
tivos 


1 — Conhecimento da origem dos sedimentos 
do estuário e da costa adjacente, suas tra- 
jectórias e respectivas condições de de- 
posição. 


2 — Conhecimento da evolução do estuário e 
da linha de costa. 
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3 — Fornecimento de elementos para a reali- 
zação, no L.N.E.C., dum modelo de fundo 
móvel do estuário do lejo, trabalho solici- 
tado pela Administração-Geral do Porto 
de Lisboa. 


4 — Conhecimento petrográfico dos sedimentos. 


5 — Contribuição para o conhecimento da Geo- 
-história do estuário do Tejo. 


Dos objectivos enunciados, os três primeiros 
são os de maior aplicação prática. Os restantes, 
embora mais académicos, não poderiam deixar 
de estar no espírito do autor, em virtude de 
constituirem contribuição para o conhecimento 
geológico de uma região importante. Com esse 
fim estudaram-se cerca de 300 amostras distri- 
buídas pelo estuário, pelas praias da margem 
direita e da margem esquerda do rio, pelas praias 
da costa para norte e para sul do estuário e por 
vários afloramentos da margem esquerda. 

O trabalho de campo decorreu em duas cam- 
panhas. A primeira, no Verão de 1962, durante 
a qual foram colhidas amostras de praia desde o 
Guincho até Algés e desde a Trafaria até cerca 
de 5km. para sul da Lagoa de Albufeira. A se- 
gunda, na Primavera e Verão de 1964, durante 
a qual foram colhidas algumas amostras da ca- 
mada superficial do leito do rio, amostras de 
afloramentos miocénicos, pliocénicos e pleisto- 
cénicos, amostras de praia na costa atlântica, 
especialmente para norte do Guincho, e amostras 
nos cachopos do norte e do sul. 

A maior parte das amostras do leito do rio, 
com excepção das colhidas na Primavera de 1964, 
e já referidas, foram postas à disposição do 
L.N.E.C. pela Administração-Geral do Porto de 
Lisboa. 

Os dados topográficos actuais relativos ao es- 
tuário e grande parte dos valores da salinidade 
foram fornecidos pela Missão do Continente e 
Ilhas Adjacentes do Instituto Hidrográfico da 
Marinha. 


2 — CONSIDERAÇÕES GERAIS 


O rio Tejo nasce em Espanha, na Serra de 
Albarracim, a 1593m de altitude e vem desa- 
guar no (Oceano Atlântico, junto a Lisboa. Tem 
uma extensão de cerca de 1100km, dos quais 
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apenas 230 em território português, e a área da 
sua bacia hidrográfica é de cerca de 86 000 km. 

Considera-se neste trabalho que o seu estuário 
(Fig. 1), se estende desde um alinhamento que 
passa junto da Póvoa de Santa Iria até à foz 
(alinhamento pelas torres de S. Julião e do Bugio), 
com área de çerca de 24000ha (Arantes e Oli- 
veira, 1941); recebe directamente águas drena- 
das por alguns afluentes, dos quais, os mais im- 
portantes são, na margem esquerda, o rio Sor- 
raia, e, na margem direita, o rio Trancão, a ri- 
beira do Jamor e a ribeira de Barcarena. A 
comunicação com o Oceano é feita por três 
barras, das quais, a mais importante é a Barra 
Grande compreendida entre dois cachopos, o do 
Norte e o do Sul, com a profundidade mínima 
de 11 m. 

Os aspectos da costa adjacente ao norte e ao 
sul do estuário, devido principalmente à diferen- 
ciação litológica das formações que a constituem, 
são manifestamente diferentes. Para sul, a costa 
é baixa e arenosa, com uma praia contínua com 
cerca de 30 km de extensão que vai terminar já 
nas proximidades do Cabo Espichel. Para norte, 
a costa é mais abrupta e recortada, com grande 
número de pequenas praias entre formações ro- 
chosas que penetram no mar. 

O estuário não se apresenta com as mesmas 
características ao longo de toda a sua extensão. 
Para tal concorrem vários factores como, por 
exemplo, a distribuição dos fundos, a salinidade, 
a geologia e a morfologia das formações que 
marginam o estuário, as marés, as correntes e a 
evolução dos fundos e da linha de costa. Embora 
o estudo experimental destes factores não tenha 
sido objecto do presente trabalho, parece indis- 
pensável apresentar algumas considerações a seu 
respeito, baseadas em elementos obtidos de várias 
fontes. 


2.1 — Distribuição dos fundos 


Da observacão dos planos hidrográficos da 
Barra e do Porto de Lisboa, verifica-se uma de- 
sigual distribuição dos fundos ao longo do es- 
tuário. Deste ponto de vista pode dividir-se o 
estuário em duas grandes zonas (Fig. 2). 

A primeira, que se inicia junto à foz e se pro- 
longa para montante do alinhamento Terreiro do 
Paço-Cacilhas, com profundidades da ordem dos 
30m. e atingindo, em certos locais, valores su- 
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Fig. 1 — Estuário do Tejo 


periores a 40m. Nesta zona, as profurdidades 
do leito do rio junto à margem esquerda são 
bastante superiores às correspondentes na mar- 
gem direita. A segunda zona, para montante da 
primeira, com profundidades inferiores a 10m, 
apresenta-se sulcada por alguns canais (Fig. 2). 
Junto das margens, em especial na margem es- 
querda, na baixa-mar o leito do rio fica a des- 
coberto. 


2.2 — Evolução dos fundos do estuário e da 
linha da costa. 


Datam do século XVI as primeiras cartas do 
estuário do Tejo, com bastante pormenor, in- 
cluindo já a indicação de profundidades (Lemos, 
1964). 

No entanto, só em 1842 se iniciou o levanta- 
mento de planos hidrográficos da barra e es- 
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tuário do Tejo que, pela sua orientação, per- 
mitem um estudo da evolução dos tundos e da 
linha de costa ao sul da barra por comparação 
com os planos mais recentes. Os planos hidro- 
gráficos mais recentes que se possuem são, para 
a Barra do Porto de Lisboa, o levantado em 1929 
e actualizado em 1960 e, para o Porto de Lisboa, 
o que se encontra actualmente em vias de con- 
clusão. 

Da análise e comparação destes planos hidro- 
gráficos, podem destacar-se, como conclusões 
importantes, as seguintes: 


a) Em relação à evolução dos fundos do es- 
tuário pode considerar-se que na sua grande 
parte (desde aproximadamente Sacavém, até um 
pouco para montante do alinhamento Belém-Ba- 
nática) tem havido pequena evolução. Na fig. 3 
apresentam-se dois perfis com o mesmo alinha- 
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Fig. 2 — Batimetria do estuário 


mento e obtidos dos planos hidrográficos de 1878 
e 1965; por ela se pode verificar que em toda a 
zona referida, no decorrer destes anos se veri- 
ficou uma sedimentação ordem de 1 a 2 metros. 

A estudar-se a evolução dos fundos, nesta 
zona, há a ter em consideração os inúmeros tra- 
balhos de dragagem que se têm vindo a rea- 
lizar. 

No início do estuário, sobretudo em certas 
zonas do perfil transversal, verifica-se, contudo, 
acentuada sedimentação de material silto-argiloso; 
essas zonas correspondem ao limite jusante dos 
mouchões. 

Nas proximidades da foz verifica-se que no 
decorrer do último século houve importante de- 
posição de sedimentos. Com efeito, observam-se 
diferenças de cota que chegam a atingir 14 m. 
Na fig. 4 apresentam-se dois perfis com o mesmo 
alinhamento obtidos dos planos hidrográficos 
de 1879 e 1960; a sedimentação, como é evi- 
dente na figura, teve lugar em quase todo o 
perfil transversal, 
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Mais adiante procurar-se-á justificar a pre- 
sença deste material nesta zona a definir a sua 
origem. 


b) Em relação à evolução da linha de costa 
atlântica, há a distinguir a costa ao norte e a 
costa ao sul do estuário. 

Ao norte, dada a litologia e a morfologia da 
costa, não é sensível, por comparação dos planos 
hidrográficos atrás mencionados, qualquer evolu- 
ção; contudo, sabe-se que as praias da baia de 
Cascais, em especial a praia do Estoril, vinham 
a sofrer nos últimos anos, emagrecimento pro- 
gressivo. A realização de um esporão, nesta 
última, no extremo oriental da praia, dificultou a 
continuidade do transporte dos sedimentos. Ao 
sul do estuário, a linha de costa tem sofrido 
grandes alterações, em especial nas proximidades 
da foz. Sobretudo nos últimos anos, tem-se assis- 
tido a um emagrecimento contínuo das praias 
para norte da Fonte da Telha e, em especial, 
pata norte da Costa da Caparica. 
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Fig. 4 — Evolução dos fundos do estuário segundo o perfil Algés-Trafaria 
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2.3 — Salinidade 


O estuário do Tejo faz parte da zona do rio 
de características marítimas, que se estende, se- 
gundo Lemos (1964), até Vila Franca de Xira. 
A sua alimentação é mista de água doce e água 
salgada. A proporção de cada uma delas varia 
com a época do ano, e pode ser avaliada a par- 
tir da salinidade média das águas do mar nas 
proximidades da barra. 

Infelizmente, as informações relativas aos va- 
lores da salinidade no estuário do Tejo são muito 
reduzidas. Serviram de base às considerações que 
se seguem os valores apresentados no trabalho 
«Esgotos de Lisboa» (Arantes e Oliveira, 1941) 
e alguns valores, ainda não publicados, forne- 
cidos pelo Instituto Hidrográfico da Marinha e 
e obtidos em Abril e Agosto de 1964, 

No quadro 1 apresentam-se os valores das sali- 
dades médias em vários pontos do estuário 
observados à superfície, no fim da estiagem. 


QUADRO 1 
| Valores de salinidade g/] 
Local ERC SR Sd 

| Minimo | Máximo | Médio 
Caxias 34,60 | 35,8 | 35,0 
Cova do Vapor 35,0 | 35,8 | 35,2 
Pedrouços 34,7 | 35,2 | 35,0 
Belém 34,2 | 35,5 | 34,9 
Alcântara 32,5 | 84,8 33,4 
Olho de Boi 33,9 | 35,0 | 34,5 
Cais do Sodré 33,2 | 34,5 | 33,9 
Margueira | 221 | 33,1 | 83,2 
Este de Cacilhas 91,6 | 35,0 | 33,8 
Montijo 91,6 | 34,8 | 33,2 
Xabregas | 31,6 | 852 | 884 
Espigão do Montijo 20,5 | 29,5 | 25,0 
Cabo Ruivo 13,0 | 3240 | 22,5 
Olivais 15,0 | 31,7 | 23,4 
Moscavide | 15,0 | 29,3 | 22,9 
Pancas | 16,7 | 237 | 20% 
Sacavém 5,0 | 28,9 18,5 
Póvoa St.“ Iria a as 3.9 


Em face destes valores médios da salinidade à 
superfície, pode dividir-se o estuário do Tejo em 
várias zonas (Fig. 5): A primeira, que se estende 
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desde a barra até às proximidades de Belém, 
com valores de salinidade próximos dos da água 
do mar; a segunda, para montante de Belém, 
que vai aproximadamente até ao alinhamento 
Terreiro do Paço-Cacilhas, com salinidades mé- 
dias da ordem de 34 g/l; a terceira, desde este 
alinhamento até um pouco para montante do 
alinhamento Xabregas-Montijo, com salinidades 
médias de 33,3 g/l; a quarta, para montante 
deste alinhamento até às proximidades de Saça- 
vém, com valores médios da ordem dos 24 g/1. 
Para montante, conhecem-se apenas os valores 
da salinidade em frente a Sacavém, da ordem 
dos 18,5 g/l, e no início do estuário, junto à 
Póvoa de Santa Iria, da ordem dos 3,3 g/1. 

Na época das chuvas todos os valores atrás 
referidos diminuem considerâvelmente, em espe- 
cial próximo da superfície, em virtude da grande 
quantidade de água doce que aflui ao estuário. 

Para montante de Alhandra, os valores mesmo 
na estiagem (Agosto 1964) são nulos. 

Pelo quadro 2, transcrito do trabalho referido 
(Arantes e Oliveira, 1941), podem verificar-se, 
para medidas obtidas a várias profundidades, as 
diferenças de salinidade em g/I, para a mesma 
estação onde se obtiveram valores na vertical, 
em dois períodos do ano climáticamente dife- 
rentes (princípios de Agosto e fins de Dezembro 
de 1937). Essas diferenças deixam quase de se 
sentir, abaixo de uma determinada profundidade, 
o que evidencia incompleta mistura das águas 
doce e salgada devido à diferença de densidades. 
Contribuindo para uma maior diferença nos va- 
lores obtidos mas duas épocas, há que ter em 
conta a evaporação que se faz sentir intensa- 
mente na estiagem. 


2.4 — Transporte dos sedimentos 
2.4.1 — Generalidades 


O transporte de sedimentos pela água, prin- 
cipal agente transportador dos sedimentos do 
estuário do Tejo e da costa atlântica, é função, 
sobretudo, da velocidade, turbulência e salini- 
dade da água e da forma, dimensão e peso espe- 
cífico das partículas transportadas. 

Habitualmente, as partículas de maior diá- 
metro são transportadas por arrastamento pelo 
fundo, as de diâmetro intermédio transportadas 
aos saltos e as de menor dimensão transportadas 
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QUADRO 2 


Profundi- Cruz Quebrada | Caxias 
dade - ' nego 
Ag. Dez. j Ag 
(m) 3 27 3 
35,46 | 22,61 
5 3544 | 30,85 25,46 
10 38,43 32,57 35,48 
20 ca 34,45 a 
25 35,48 34,47 a 
Fundo 30,41 34,67 35,57 


28,68 


35,46 22.927 35,46 | 26,22 | 35,41 | 26,08 
32,57 


o 


Ag. Dez. 


35,50 | 82,12 | 85,50 | 27,9 

35,52 | 34,29 | 85,55 | 3481 

34,67 sa 35,85 

cos ge si sa 35,89 
34,76 35,57 | 35,80 | 35,64 | 35,44 


em suspensão. Contudo, as variações de velo- 
cidade e de turbulência podem sujeitar, em re- 
gimes diferentes, a mesma partícula aos três 
tipos de movimento mencionados. 

No caso do estuário do Tejo, interessa consi- 
derar a acção transportadora do rio e a acção 
transportadora do mar. 


2.4.2 — Acção transportadora do rio 


A capacidade de transporte fluvial varia con- 
siderávelmente com as épocas do ano, sendo 
função do caudal do rio. Na mesma época varia, 
contudo, ao longo do perfil longitudinal, dimi- 
nuindo com a aproximação da foz, em especial 
se, como no caso do Tejo, o rio corre na sua 
parte final em regiões mais ou menos planas. 

Num estuário, além das correntes originadas 
pelo caudal fluvial há a considerar as correntes 
de maré. No estuário do Tejo essas correntes 
são, com certeza, muito intensas, como deixam 
prever os valores das amplitudes de maré que se 
conhecem; 4,5m de amplitude máxima junto ao 
Montijo, 3,75 m junto a Paço d'Arcos e osci- 
lando pelos 4,3m na maior parte do estuário. 

A água salgada, mais densa, tem tendência a 
penetrar no estuário pelo fundo, nunca se mis- 
turando completamente com a água doce. Este 
facto ocasiona correntes de densidade que pro- 
vocam turbulências que impedem, dentro de 
certa medida, a deposição dos materiais trans- 
portados pelo rio em suspensão e que vão sendo, 
assim, arrastados até ao mar. Também a inversão 
da corrente de maré, por não se verificar simul- 
tâneamente em todo o perfil vertical, contribui 
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para evitar a deposição de sedimentos em quase 
todo o estuário. Com efeito, medições de velo. 
cidades, próximo da margem esquerda, na secção 
do vale de Alcântara, a montante de Belém, nos 
nos meses de Setembro, Outubro, Novembro e 
Dezembro de 1956, conduziram na zona referida 
às seguintes conclusões (Abecasis e Lencas- 
tre, 1957): 


a) As velocidades da corrente a 2m do fundo 
são sempre um pouco maiores na enchente do 
que na vazante (para um coeficiente de maré 70, 
na vazante a Velocidade é cerca de 70cm/s e 
na enchente cerca de 90 cm/s; 

b) As velocidades da corrente, na enchente 
e na vazante, a 2m da superfície, têm valores 
aproximadamente iguais. 

Durante a época das chuvas, em que o caudal 
sólido aumenta consideravelmente, pode obser- 
var-se, no mar, junto à foz do estuário, uma man- 
cha de cor amarela que sai a barra e penetra 
alguns quilómetros no oceano. O caudal sólido 
transportado, consequência da erosão continental 
é pois arrastado para o largo, até que aí encon- 
adequadas condições de deposição nos fundos 
marinhos, onde já não se fazem sentir nem as 
acções das correntes nem das ondas. 

Junto às margens do estuário, em especial no 
corredor entre Cacilhas e Trafaria, na margem 
esquerda, o rio vai erodindo as formações que as 
constituem, transportando, embora a pequena 
distância, o material destacado. 

Parece poder concluir-se que, no estuário do 
Tejo, actualmente, o transporte apreciável de sedi- 
mentos pelo rio se faz apenas em suspensão e 
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Fig. 5 — Distribuição de salinidade no estuário 


a maior parte desse material se vai depositar no 
mar, em zonas afastadas da foz. O material are- 
noso transportado pelo rio, por arrastamento, 
sedimenta, na sua quase totalidade, a montante 
da zona que se considerou como estuário. 

O material arenoso resultante da erosão da 
margem esquerda sofre, como se verá em 6.2.1, 
apenas, localmente, pequeno transporte. 


2.4.3 — Acção transportadora do mar 


O transporte marinho é devido, fundamental- 
mente, à acção das correntes marinhas e das on- 
das. 

Das correntes marinhas, devem ser as corren- 
tes de maré aquelas cuja intensidade mais se faz 
sentir junto à costa. O regime de correntes (su- 
perficiais e profundas) na nossa costa é contudo, 
quase desconhecido. Num esquema muito geral 
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(Freire de Andrade, 1937) pode considerar-se a 
costa Portuguesa sujeita à acção de duas corren- 
tes superficiais relativamente pouco intensas e 
possivelmente derivadas da corrente do Golfo; 
a primeira, de direcção aproximadamente nor- 
te-sul, correndo ao longo da costa norte até 
as proximidades de Peniche; e a segunda, 
de direcção aproximadamente sul-norte, e sen- 
tir-se-ia na costa, a sul de Peniche. Ao en- 
contrarem-se, as duas correntes inflectiriam para 
oeste Uma outra hipótese, mais aceitável, tam- 
bém citada por Freire de Andrade (1937) e por 
Kuenen (1950), refere-se à existência de uma 
corrente superficial, de direcção geral aproxima- 
damente N-S, derivada da corrente do Golfo, 
infiectindo para leste junto da ponta de Sagres e 
seguindo paralelamente à costa algarvia e à costa 
de Espanha até Gibraltar. Nas proximidades da 
foz do Tejo, devido a movimentos turbilhonares 
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resultantes da morfologia da costa, essa corrente 
ramificar-se-ia para leste, indo incidir na costa 
atlântica, ao sul do estuário, com direcção apro- 
ximadamente sudoeste-nordeste. 

Estas correntes superficiais não exercem, con- 
tudo, acção apreciável no transporte dos sedi- 
mentos. 

Quanto às correntes de fundo, nas proximida- 
des da foz do estuário, não se sabe até que ponto 
exercerão acção de transporte sobre os sedimen- 
tos, pois, com os métodos utilizados neste 
estudo, apenas foi possível considerar a re- 
sultante da acção conjunta das correntes e das 
ondas. 

As ondas são dum modo geral responsáveis 
pelo transporte dos sedimentos para as praias e, 
consequentemente, responsáveis pela sua evolu- 
ção. São só os sedimentos existentes entre a 
praia e uma dada cota, que pode variar com as 
características locais, que estão em constante mo- 
vimento de vai-vem, numa direcção correspon- 
dente à direcção de propagação da onda. Os res- 
tantes, a não se verificar a acção de correntes 
profundas, permanecerão imóveis até que se mo- 
difigquem as condições. Se as características das 
ondas não se modificarem apreciâvelmente, a sua 
acção de vai-vem sobre esses sedimentos resulta 
numa boa calibração das partículas que os cons- 
tituem, pois evita que, nessa zona, sedimente 
material de diâmetro inferior a uma determinada 
dimensão ; por outro lado, a velocidae caracteris- 
tica desse movimento é inferior à necessária para 
transportar, para o local, material de dimensão 
superior. 

Sempre que a direcção de propagação da onda 
é mais ou menos oblíqua à linha da costa, e isso 
é o caso mais frequente, resulta daí sobrepôr-se 
ao movimento de vai-vem na direcção da onda, 
um movimento longitudinal dos sedimentos ao 
longo da costa. 

Verifica-se (ver por exemplo os trabalhos da 
Comissão de Estudo de Erosão Costeira ao sul 
do Tejo) que nas proximidades da foz do estuá- 
rio, na costa atlântica sul, a direcção de propa- 
gação das ondas, nas imediações das praias, é 
aproximadamente sudoeste-nordeste, em especial 
na zona entre a Caparica e a Cova do Vapor, o 
que ocasiona forte transporte longitudinal de 
material arrancado às respectivas praias. 

Nas praias da margem norte do estuário, e 
desde Cascais até S. Julião da Barra, observa-se 
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também, transporte longitudinal, mas de ocidente 
para oriente. 


2.5 — Geologia 


A fig. 6 representa um esboço geológico das 
formações que marginam o estuário e a costa 
atlântica adjacente, feito com base nas folhas 
n.º 1,2, 3 e 4 da Carta Geológica dos Arredo- 
res de Lisboa na escala 1:50 000 e na Carta Geo- 
lógica de Portugal na escala 1:500 000. 

Todas as formações são de natureza sedi- 
mentar, com excepção da formação basáltica 
que aflora na margem direita, aproximadamente 
entre Lisboa e Paço d'Arcos, e das formações 
pertencentes à intrusão da Serra de Sintra. 

Na margem direita do estuário afloram, de 
montante para jusante, formações do Miocénico 
marinho, constituídas por areias, siltes e argilas 
que contactam com as rochas da formação basál- 
tica de Lisboa a jusante do Terreiro do Paço. 
Desde o início do estuário até um pouco a ju- 
sante do rio Trancão, o Miocénico encontra-se, 
nas proximidades da margem, coberto por for- 
mações recentes. Os basaltos e tufos vulcânicos 
da formação basáltica estendem-se, praticamente 
sem interrupção, até Paço d'Arcos, onde con- 
tactam com os calcários do Cretácico, que se 
vêem na costa ao norte do estuário até um pouco 
para além da praia grande do Guincho. Aí, o 
Cretácio contacta com os calcários e «xistos» do 
Jurássico que afloram até ao contacto com as 
rochas plutónicas da Serra de Sintra. Estas aflo- 
ram até um pouco para norte do Cabo da Roca, 
onde contactam, de novo, com as formações do 
Jurássico. Para norte da Praia da Adraga, e prá- 
ticamente sem interrupção até às proximidades 
de Ribamar, afloram novamente as formações do 
Jurássico. 

Na margem esquerda do estuário, a montante 
de Cacilhas, predominam as formações do Pliocé- 
nico cobertas, em grande extensão, por formações 
recentes. Nas proximidades de Cacilhas e até à 
Trafaria, próximo da foz, constituindo a margem 
esquerda do canal, aflora uma escarpa miocénica 
sulcada por alguns vales pouco largos. Os es- 
tratos que a constituem têm inclinação geral para 


sul da ordem dos 6 ou 7º (Oliveira, 1964), o que 
dentro de certa medida, contribui para a relativa 


estabilidade da escarpa. 
Esssa inclinação, que no perfil de implantação 
da Ponte sobre o Tejo é de 7º para SSW, corres- 
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ponde aproximadamente à inclinação do «bed- 
-rock” (formação basáltica) (Moitinho de Al- 
meida, 1962). Devido à ausência de fenómenos 
tectónicos importantes, essa inclinação mantem-se 
mais ou menos constante em toda a escarpa. 
Desde a Trafaria até à foz, marginam o es- 
tuário formações dunares recentes. Toda a costa 


Sezimbra e um outro do Miocénico nas proxi- 
midades da Arrábida. 

O estudo geológico para a implantação da 
Ponte sobre o Tejo evidenciou que o vale, na 
parte terminal do rio, não é de origem tectónica, 
como alguns autores pensaram, mas sim um vale 
de erosão. No perfil da Ponte, o “bed-rock” en- 
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FORMAÇÃO BASÁLIICA 


ROCHAS PLUTONICAS 
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Fig. 6 — Esquema geológico das formações que marginam o estuário e a costa atlântica adjacente 


atlântica ao sul do estuário, desde a sua foz 
até um pouco para o sul da Lagoa de Albu- 
feira, é constituída por formações recentes 
que contactam, para o interior, com o Miocé- 
nico e o Pliocénico. Para sul afloram, junto 
à costa, as rochas do Miocénico que contactam 
com as formações do Cretácico. Nas proximi- 
dades do Cabo Espichel afloram as rochas do 
Jurássico, que se estendem, para oriente, até às 
proximidades de Setúbal interrompidas por um 
pequeno afloramento da formação basáltica em 
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contra-se, nas proximidades da margem sul, à 
cota — 80 m, subindo progressivamente para norte 
até — 25m, voltando a descer ligeiramente nas 
proximidades da margem norte. Segundo Zu- 
zarte Mendonça, citado por Freire Andrade (1937), 
as profundidades máximas do “bed-rock” no 
perfil Beato-Montijo são da ordem dos 70m, 
Estes números conduzem à conclusão que o nível 
médio das águas do mar estava aproxidamente 
uma centena de metros abaixo do actual, quando 
terminou a escavação da zona terminal do vale 
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PLANTA 1 - LOCAIS DE COLHEITA DAS AMOSTRAS 


10 km 


do Tejo, no final da grande regressão do Pleis- 
tocénico. 


3 —- MÉTODOS UTIZADOS 
3.1 — Descrição dos métodos 
3.1 1 — Colheita de amostras 


As amostras de praia foram colhidas em perfis 
perpendiculares à linha de costa. Em cada perfil 
colheram-se três amostras, aproximadamente nas 
linhas de baixa-mar, de preia-mar e na zona in- 
termédia. 

Em cada uma das praias situadas na margem 
norte do estuário e na costa norte, conforme a 
sua extensão, foram localizados dois ou três perfis; 
nas praias da margem sul do estuário e da costa 
sul, localizaram-se 24 perfis, sendo 12 desde a 
Trafaria até à Costa da Caparica, aproximada- 
mente equidistantes e os restantes 12 para sul 
da Costa da Caparica nas zonas de mais fácil 
acesso. 

Todas as amostras de praia foram colhidas 
com um amostrador cilindrico de cobre. Nos 
locais de colheita introduziu-se o amostrador até 
à profundidade de 30 cm pela queda de um peso 
de 10kg sobre o amostrador; obteve-se assim, 
para cada local, cerca de 1 kg de amostra. 

Quase todas as amostras do leito do rio foram 
retiradas de tarolos de sondagens realizadas, em 
várias épocas, pela A. G.P.L. Além das amostras 
relativas à camada superior do leito do rio, foram 
retiradas mais uma ou duas de cada tarolo, de 
camadas mais profundas. Procurou-se serem estas 
amostras tão representativas quanto possível, 
homogeneisando todo o material existente de 
cada camada, retirando em seguida uma parte 
que se julgou suficiente para os ensaios a rea- 
lizar. 

As restantes amostras do leito do rio, bem 
como a dos cachopos, foram retiradas com um 
amostrador do tipo Van Veen. 

As amostras de afloramentos miocénicos, plio- 
cénicos e pleistocénicos da margem esquerda 
foram colhidas com uma colher, em cortes na- 
turais previamente limpos, tendo-se procurado 
que elas fossem representativas da formação. 

Na Planta 1 estão indicados os locais de co- 
lheita de amostras estudadas do ponto de vista 
mineralógico. 


3.1.2 — Minerais pesados 


O método usado no estudo do minerais pe- 
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sados (d >> 2,89) foi o desenvolvido pela escola 
holandesa e descrito por Doeglas (1940); indi- 
cam-se em seguida, resumidamente, as opera- 
ções realizadas. Separou-se uma porção repre- 
sentiva de cada amostra, que se tratou a quente 
com HCl concentrado e, em seguida, com HNO; 
concentrado ; deste modo eliminaram-se os car- 
bonatos, os sulfuretos e a matéria orgânica e 
clarificaram-se os grãos para serem estudados 
ao microscópio. Depois de tratada, cada amostra 
foi passada nos peneiros de malha 0,50 e 0,05mm. 
O estudo dos minerais pesados foi realizado ape- 
nas na fracção 500-50 4, e a sua separação feita 
em bromofórmio (d==2,89). Os minerais pe- 
sados, assim obtidos, foram montados em bál- 
samo do Canadá para serem identificados e con- 
tados ao microscópio polarizante. Foram identi- 
ficados em cada lâmina cerca de 100 grãos, não 
opacos e não alterados, e em seguida calculada 
a sua percentagem. A percentagem dos opacos, 
alterados e fragmentos de rochas foi calculada 
em relação ao número total de grãos contados. 


3.1.3 — Granulometria 


Para determinar a granulometria das amostras 
usaram-se os métodos de peneiração e pipetagem, 
respectivamente para as fracções arenosa e siltosa. 

A granulometria das amostras de praia foi 
obtida apenas por peneiração, visto a fracção 
silto-argilosa ser, praticamente, inexistente. 

A granulometria das amostras do rio foi obtida 
pela combinação dos dois métodos referidos. 

A descrição pormenorizada dos métodos está 
feita por vários autores, pelo que apenas se indi- 
cará, resumidamente, como se procedeu. 


— Granulometria por peneiração 


As amostras de praia foram lavadas com água 
destilada, para eliminar o NaCl e posterior- 
mente secas numa estufa. A totalidade de cada 
amostra seca foi peneirada manualmente nos 
peneiros n.º 5 (4,00 mm) e 8 (2,38 mm) da 
série A.S.T.M. A fracção passada no peneiro 8 
foi então dividida por esquartelamentos sucessivos 
até se obterem entre 50 e 100 gramas de amostra 
representativa. Esta porção foi, então, peneirada 
durante 15 minutos numa máquina Endecott's 
nos peneiros A.5.T.M. nº” 14, 20, 35, 50, 80, 
140 e 230 de malhas respectivamente 1,41 mm, 
0,84 mm, 0,50 mm, 0,297 mm, 0,177 mm e 
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e 0,062 mm, sendo as fracções retidas nesses pe- 
neiros devidamente pesadas. 


— Granulometria por peneiração e pipetagem 


Cada amostra foi desagregada e em seguida 
peneirada manualmente no peneiro A.S.T.M. 
n.º 10 (0,787 mm). Foram pesados çerca de 25 
gramas da fracção passada no peneiro n.º 10, 
devidamente homogeneizada, aos quais se juntou 
água oxigenada até completo ataque da matéria 
orgânica, e em seguida HCl para serem remo- 
vidos os carbonatos; cada amostra foi então 
lavada num funil de Biichner, juntando-se de- 
pois 50c.c. de solução de oxalato de sódio (dis- 
persante); em seguida foi lavada num peneiro 
de malha 0,050 mm tendo a fracção que passou 
no peneiro sido recebida numa proveta gra- 
duada, à qual se juntou água destilada até per- 
fazer 1000c.c. O resíduo que ficou nesse pe- 
neiro foi seco e peneirado mecânicamente, durante 
15 minutos, nos peneiros A.5.T.M. n.º 14, 20, 
35, 50, 80, 140 e 230. Para determinação da 
granulometria da fracção < 0,050 mm, em pri- 
meiro lugar, agitou-se a suspensão para a homo- 
geneizar e retiraram-se imediatamente 20c-.c. 
com a pipeta. O peso total da suspensão foi cal- 
culado a partir do peso dessa fracção inicial- 
colhida. Juntando o peso da suspensão com o 
da fracção >> 0,050mm obteve-se o peso total 
do sedimento, em relação ao qual se fez o cál- 
culo das percentagens das diferentes fracções 
granulométricas. 

Determinaram-se em seguida os pesos das 
fracções < 0,032, < 0,016, < 0,008, < 0,004 e 
<. 0,002 mm por pipetagens sucessivas, aumen- 
tando progressivamente o intervalo entre o agi- 
tar da amostra e a colheita dos 20 c.c. da sus- 
pensão. 

O método da pipeta baseia-se na lei de Stokes, 
que permite calcular o diâmetro das partículas 
em função da sua velocidade de sedimentação. 


y — 2 (di — da) gr? 
97mn 


(1) 


sendo 


v — velocidade da partícula em cm/s 
di — densidade da partícula 
ds — densidade do fluído 


g — aceleração da gravidade em cm/s? 
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r —raio da partícula em cm 
n — viscosidade do fluído 


Sendo conhecidos os valores da densidade das 
partículas, da densidade e da viscosidade da 
água e da aceleração da gravidade, a fórmula (1) 
pode representar-se: 


v=Cr' (2) 
em que 
C — 2(di — dog 
9 n 


Os valores de C estão calculados em tabelas. 
3.1.4 — Carbonatos 


Para a determinação da percentagem dos car- 
bonatos aplicou-se o método de Scheibler, que 
consiste na medição do volume de CO» resul- 
tante da reacção de HCI concentrado com uma 
pequena porção de amostra desagregada, após 
cerca de duas horas de agitação mecânica. 

Esse volume é comparado com o volume de 
CO: libertado na reacção da mesma quantidade 
de HCl com uma solução de NasCO; (padrão). 
O resultado é expresso em percentagem de car- 
bonatos. 


3.2 — Precisão no doseamento dos minerais 
pesados 


No estudo dos minerais pesados dos sedimen- 
tos há que considerar os erros cometidos na 
colheita das amostras e os erros laboratoriais. 
Estes compreendem os erros cometidos na sepa- 
ração por líquidos pesados, na montagem em 
lâmina e no doseamento dos minerais pesados. 

Os autores são, dum modo geral, unânimes 
em atribuir a maior percentagem de erro labora- 
oriasao doseamento dos minerais pesados, dando 
especial relevo ao número de minerais pesados 
contados em cada amostra. 

A maioria dos autores mais recentes usa con- 
tar em cada lâmina cerca de 100 grãos não opa- 
cos, seguindo o método descrito por Doeglas 
(1940). No entanto, muitos outros estimam nú- 
meros muito superiores para que a precisão seja 
satisfatória. Assim, H. Deeke conta 300 grãos, 
V. Leinz 400 a 600. O. Voelcker 500, K. H. Sin- 
dowsky 400, W. Scheidbarer 1 000. Estes autores 
são citados por Van Andel (1950), que conta 
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100 grãos não opacos. Duplaix (1958) conta to- 
dos os grãos montados (cerca de 200 em cada 
lâmina). 

É evidente que a precisão aumentará desde 
que se aumente consideravelmente o número de 
grãos contados. No entanto, o duplicar ou tripli- 
car esse número equivale praticamente a dupli- 
car ou triplicar o tempo que demora a realiza- 
ção do trabalho, sendo os resultados obtidos ape- 
nas um pouco menos dispersos e uma das razões 
que leva o estudo dos minerais pesados a ser 
empregue por grande número de autores é, pre- 
cisamente, o facto de se obterem rapidamente 
resultados satisfatórios. 

No presente trabalho foram contados, para 
cada lâmina, cerca de 100 grãos não opacos e 
não alterados, tendo-se calculado as percenta- 
gens dos grãos opacos e dos alterados em rela- 
ção ao número total de grãos contados. Os alte- 
rados foram contados separadamente por atingi- 
rem, por vezes, percentagens muito elevadas, o 
que, a incluí-los nos 100 grãos não opacos, iria 
afectar bastante o valor das frequências dos ou- 
tros minerais. 

No quadro 3 figuram as percentagens dos dife- 
rentes minerais pesados das amostras n.º5 80, 
211 e 296, contando para cada uma e na mesma 
lâmina 100, 200, 300, 400 e 500 grãos. Na 
última coluna apresentam-se as diferenças entre 
as percentagens de cada mineral entre a conta- 
gem de 100 e a contagem de 500 grãos. Pode 
verificar-se que as diferenças máximas observa- 
das são da ordem dos 5 “/o, valor que é inferior 
ao que normalmente se exige para que se possa 
considerar um mineral como característico de de- 
terminada associação. Além disso, as conclusões 
baseiam-se habitualmente, em diferenças muito 
superiores. 

Outros aspectos a considerar, quando se estu- 
dam os minerais pesados de uma dada região, 
dizem respeito aos limites granulométricos entre 
os quais se deve fazer o estudo dos minerais pe- 
sados e ao número de fracções granulométricas 
a estudar em cada amostra. 

A fixação desses limites deve ser feita no iní- 
cio de cada trabalho e para cada caso. Com 
efeito, a dimensão dos grãos que constituem o 
sedimento é fundamentalmente função da dimen- 
são dos minerais na rocha mãe e das condições 
de transporte; no entanto, dado que, em condi- 
ções normais, sedimentam ao mesmo tempo as 
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particulas de igual equivalente hidráulico, não é 
de prever o aparecimento de minerais pesados 
acima de determinadas dimensões. As opiniões 
variam, no entanto quanto à dimensão limite do 
estudo. Parece, dum modo geral, satisfatório o 
limite superior de 500 |, mais habitualmente 
usado, pois acima deste valor a frequência dos 
minerais pesados é muito reduzida; por outro 
lado, não é de esperar o aparecimento de mine- 
rais pesados na fracção argilosa, e mesmo na 
siltosa a sua frequência é habitualmente muito 
reduzida. A maior parte dos autores limita infe- 
riormente em 50 ou 60 a fracção a estudar. 
Com efeito, abaixo destes valores o estudo mi- 
croscópico e até a separação no bromofórmio 
começam a revestir-se de maior dificuldade e 
resultam num muito maior dispêndio de tempo. 
Há autores, e especialmente os que se dedicam 
à mineralogia de solos, que levam o estudo até à 
fracção 20 ». Também Van Andel (1954) estu- 
dou os minerais pesados nas fracções 75-50 pg, 
50-30 u e 30-10 » em algumas das amostras 
mais finas. Conclui, nesse trabalho, ser de 30 
o limite inferior aproximado para o aparecimento 
dos minerais pesados. Esse facto não é no entanto 
suficiente para levar o estudo até esse limite, 
pois os minerais que aparecem na fracção 50-30 p 
não são muito significativos. A maior parte da 
fracção é constituída por zircão, que é um mine- 
ral característico da fracção fina dos sedimentos, 
Com o fim de melhor definir as fracções mais 
representativas nos sedimentos do Tejo e da costa 
adjacente, separaram-se várias amostras nas frac- 
ções compreendidas entre os limites 840-500 p, 
500-297 wu, 297-177 wu, 177-105 4, 105-62 u e 
62-50 Lu. 

No quadro 4 figuram as percentagens dos mi- 
nerais pesados de cada fracção em relação ao 
peso total de minerais pesados de cada amostra. 
Nas últimas colunas figuram, respectivamente, 
os valores das medianas da amostra total (frac- 
ções leve e pesada) e dos minerais pesados. 

Conforme se pode verificar, nas fracções 
840-500 wu e 62-50 4, as percentagens dos mine- 
rais pesados são, dum modo geral, muito peque- 
nas em relação às das outras fracções. Na fracção 
840-500 » das amostras 5 e 119, as percentagens 
são elevadas devido à presença de grande quan- 
tidade de alterados e fragmentos de rochas. Em 
algumas das amostras do rio que contêm muito 
material fino, a percentagem de minerais pesados 
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QUADRO 4 


0/y de minerais pesados nas fracções 


q E 4 s. s. 

bg o a R 

Z E p q T 

o o) Es 

EA he EA 
5 40,3 33,9 19,3 
35 8,0 13,6 45,1 
59 6,0 9,0 38,8 
so 0,9 39,9 52,9 
119 44,7 31,0 10,4 
131 1,2 38,8 28,0 
158 0,0 6,5 26,1 
185 04 | 1,8 18,2 
188 0,0 0,7 1,0 
205 ii 74 48,6 


so 

es “ss 

é 3. 3. E E à 
ç ; g * | 8% 
R ê 3 Zs 
5,8 0,8 0,4 480 u 320 u 
9,4 2,5 1,2 308 u | 207 p 
40,2 5,0 1,0 278 u. 180 pu 
5,8 0,3 0,2 259 u 280 u 
3,9 T,l 3,5 970 p 320 u 
17,2 8,6 0,2 392 q 265 pg 
45,6 17,4 4,4 4 u 150 u 
43,0 27,1 10,0 25 yu 125 pg 
35,9 50,0 12,5 8u 94 u 
24,6 13,8 9,2 308 u 180 p 


na fracção 62-50 « é apreciável, pelo que se 
resolveu, no presente trabalho, fazer o estudo dos 
minerais pesados na fracção 500-50 p. 


Parece, dum modo geral, mais razoável, para 
cada amostra, fazer o estudo duma só fracção 
(no caso presente 500-50 ), em vez de dividir a 
amostra em duas ou três fracções. Duplaix (1958) 
estuda sistemáticamente as fracções 500-350 p, 
350-160 |, 160-50 » de todas as amostras, para 
que a fracção pesada seja mais homogénea, quer 
na composição quer na granulometria, o que, se- 
gundo a autora dá maior precisão e facilita a 
montagem no bálsamo de Canadá. Pureza e 
Araújo (1956) dividem todas as amostras nas 
fracções 350-177 |» e 177-62 p, que estudam 
separadamente, para que a separação dos mine- 
rais em leves e pesados seja mais correcta. 

As razões apresentadas, pelos autores citados, 
para justificarem o estudo mineralógico das 
amostras em várias fracções granulométricas, são 
de considerar. No entanto a divisão sistemática 
de todas as amostras em fracções, só seria real- 
mente vantajosa se não houvesse interesse em 
conhecer a percentagem de cada espécie mi- 
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neral na totalidade da amostra. Desde que, para 
conhecer essas percentagens, é preciso calculá- 
-las a partir das percentagens parciais de cada 
espécie mineral em cada fracção granulométrica, 
o erro que se introduz nesse cálculo é, com cer- 
teza, muito superior ao erro cometido na sepa- 
ração dos minerais leves e pesados duma amostra, 
representada por uma só fracção granulométrica. 

Além deste facto, é de ter em conta o grande 
aumento de trabalho e consequente dispêndio de 
tempo motivados pelo grande número de separa- 
ções e montagens a fazer, de lâminas a estudar 
ao microScópio e de cálculos a realizar. Parece, 
pois, vantajoso o estudo sistemático dos minerais 
pesados, na fracção 500-50 p.. 

Contudo, e devido fundamentalmente à influ- 
ência das variações granulares, se as conclusões 
se basearem apenas nos resultados da análise das 
amostras numa só fracção, podem conduzir a 
uma interpretação incorrecta. À influência dessas 
variações granulares e as medidas a tomar para 
as atenuar, ou mesmo eliminar, estão menciona- 
das em 4.4. 


(Continua) 


